Thermocapillary Flow on Superhydrophobic Surfaces 
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A liquid in Cassie-Baxter state above a structured superhydrophobic surface is ideally suited for surface 
driven transport due to its large free surface fraction in close contact to a solid. We investigate thermal 
Marangoni flow over a superhydrophobic array of fins oriented parallel or perpendicular to an applied tem- 
perature gradient. In the Stokes limit we derive an analytical expression for the bulk flow velocity above 
the surface and compare it with numerical solutions of the Navier-Stokes equation. Even for moderate 
temperature gradients comparatively large flow velocities are induced, suggesting to utilize this principle 
for microfluidic pumping. 



Introduction - Microtextured surfaces have mainly re- 
ceived attention due to their wetting properties [1]. A liq- 
uid drop placed on a suitably structured hydrophobic sur- 
face will only be in contact with the material on protruding 
tips, while gas is trapped in the valleys in between. In this 
so called Cassie-Baxter state nearly perfect hydrophobicity 
can be obtained, reflected in contact angles close to 180°. 
Recently, such surfaces have gained interest with respect to 
their ability for drag reduction 0, [3] and surface induced 
transport 0,0], in particular electroosmotic and diffusioos- 
motic flow. 

In this letter we analyze temperature induced Marangoni 
convection as a driving force for fluid transport along mi- 
crotextured surfaces. In particular, we focus on finned sur- 
faces as sketched in figure [Q with a temperature gradi- 
ent along or perpendicular to the fins, and the liquid be- 
ing in the Cassie-Baxter state. We use an integral relation 
for the Stokes equation to derive an analytical formula for 
the macroscopic flow velocity observed at some distance 
above the surface. Both situations are further investigated 
by numerical solutions of the Navier-Stokes equation and 
compared to the analytical formula. Our analysis focuses 
on substrates of high thermal conductivity, such as silicon. 

Fluid actuation and transport are core functionalities in 
many microfluidic systems. The most prominent examples 
for the corresponding driving mechanisms are pressure- 
driven and electroosmotic flow. Our analysis shows that 
moderate temperature gradients of the order of 10 K/cm 
can lead to fluid velocities of several mm/s for water based 
systems on superhydrophobic surfaces. Thermocapillary 
convection may thus add to the portfolio of actuation prin- 
ciples in microfluidic settings and may even enable larger 
flow velocities than typically achieved with electroosmo- 
sis. 

Marangoni flows - The stress on a liquid-gas interface 
due to a gradient in surface tension is ]6|] 



-tidier, 



(1) 



where n,j and ij are components of the interface normal and 
tangential vectors in i direction, and we use the convention 
to sum over repeated indices, are the components of the 
stress tensor and a is the surface tension. For an incom- 
pressible Newtonian fluid characterized by the viscosity rj 



the total stress tensor is 



ri(diUj +djUi) -pdi 



(2) 



consisting of a viscous part proportional to the shear rate 
tensor and a part originating from the pressure field p. The 
equation of motion for the fluid is the stationary Navier- 
Stokes equation 



d j (pu j u i 



dr ■ 



diUi = 0, 



(3) 



where the fluid is assumed to be incompressible. The left 
hand side of the momentum equation becomes negligible 
at small velocities (Reynolds numbers), and we will call 
this limit the Stokes limit, with the corresponding equation 
of motion being the Stokes equation. On the fin surface we 
assume a no-slip boundary condition, while eq. dU consti- 
tutes the corresponding boundary condition on the liquid- 
gas interface. 

The temperature in the fluid is governed by the stationary 
limit of the energy equation 



pCpUidiT = kAT, 



(4) 



where p, c p and k denote the density, specific heat capacity 
and thermal conductivity of the liquid, respectively. In or- 
der to tackle the problem analytically, we assume the den- 
sity, viscosity, heat capacity and thermal conductivity to be 
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FIG. 1 : (color online) Sketch of the geometry. The z-axis is al- 
ways chosen along the main temperature gradient, the y-axis al- 
ways normal to the structured surface. In this example we assume 
a temperature gradient along the fins (longitudinal). 
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independent of temperature while linearizing the temper- 
ature dependence of the surface tension |7J]. As long as 
the temperature differences do not become too large this is 
a suitable approximation. Moreover, we will restrict our 
analysis to a flat liquid-gas interface. 

Equating the viscous part of the stress tensor, ~ r)U/L, 
with the Marangoni stresses, ~ ^d z T, yields a dimen- 

sionless number, f3 = 3^1% , which for a given geometry 
characterizes the flow problem together with the Reynolds 
number, Re = pUL/r], and Prandtl number, Pr = c v r\jk. 
Below, we will demonstrate that (3 is closely related to the 
macroscopic slip length in a simple shear flow over tex- 
tured geometries. 

Longitudinal fins - Consider a geometry where the tem- 
perature gradient in the bulk of the substrate has the value 
(d z T) and is parallel to the fins, which defines the z- 
direction (figure [Q. The geometry is implied to be of in- 
finite extent in the plane of the substrate and at a height 
H above the substrate a symmetry plane is assumed. In 
that case there is translational symmetry in z-direction and 
hence the z-component of the temperature gradient will 
have the same value everywhere in the geometry. More- 
over, the velocity and pressure in the fluid will not depend 
on this coordinate. This simplifies the energy equation © 
to a 2D equation 

pc p (u ± V ± T + u z (d z T)) = kA ± T, (5) 

where we have used the shorthand notation u^V^ = 
u x d x + Uydy and Aj_ = d% + dy. The stationary Navier- 
Stokes equations partially decouple into a 2D perpendicu- 
lar equation for = (u x ,u y ) and a convection diffusion 
equation for the longitudinal velocity u z : 

p(u_lV_l)u_l = i]A ± u ± - Vj_p, 
p(u ± V ± )u z = r]A ± u z , (6) 
d x u x +d y Uy = 0. 

In the Stokes limit the left hand sides can be set equal to 
zero and the equation for u z decouples from the perpen- 
dicular parts. 

As noted above, the z-component of the temperature gra- 
dient along the fins is constant everywhere due to the trans- 
lational symmetry and this leads to a main flow along the 
fins towards the colder regions. Additionally, temperature 
variations perpendicular to the main flow occur, the tem- 
perature being lower close to a fin than in the middle be- 
tween two fins. Due to the symmetry of the problem, these 
perpendicular temperature gradients will only lead to per- 
pendicular vortices and not produce any net flow. 

Lorentz reciprocal theorem: - This integral relation for 
the Stokes equation [ 8] has in microfluidic applications al- 
ready proven useful in the case of electroosmotic flow 190 . 
It is obtained by assuming two solutions for the Stokes 
equation d3), one characterized by "unhatted" fields (u,p) 
and the other denoted by the respective "hatted" values 



(u,p) corresponding to some reference flow field, pos- 
sibly subject to different boundary conditions. Multiply- 
ing equation (|3]l in the Stokes limit by u, integrating over 
a control volume Q, integrating by parts and subtracting 
the corresponding equation where "hatted" and "unhatted" 
fields are interchanged, leaves only the boundary integrals 

i] / dA^UjdiUj-UjdiUj) = / dA^u.p-u^), (7) 
Jon Jan 

where the continuity equation was used to simplify the 
pressure term. 

Let us now take as the reference flow a Couette flow over 
the fin geometry, fulfilling the no-slip boundary condition 
at the solid-liquid and a zero shear stress boundary condi- 
tion at the liquid-gas parts of the interface. Far away from 
the substrate a constant shear rate 7 is assumed. This flow 
has been studied by Philip II 1011 . who gives an analytical 
formula for the flow field and the corresponding macro- 
scopic slip velocity, expressing that far away from the sub- 
strate it appears as if the liquid glides over the microstruc- 
tured surface. However, we shall not need the details of this 
flow field but only the "global" parameter, the slip length 
$1 over the surface. In particular, the flow field far away 
from the surface has the form 

w~ 7 (y + /3 ( )e z . (8) 

As our control volume we take a rectangular box span- 
ning the distance from the center of some fin to the center 
of its neighbor, stretching an arbitrary distance along the 
fins and being sufficiently high, such that the flow field on 
the top surface is unaffected by the details induced by the 
bottom surface, the plane containing the solid-liquid and 
liquid-gas interfaces. Assuming the top surface to be lo- 
cated at a distance H away from the bottom surface we find 
that at this position the Couette flow is approximately given 
by equation ©, while for the ("unhatted") thermocapillary 
flow a vanishing shear stress can be assumed, since the 
driving force is surface tension. This leads to a constant 
flow velocity at at the top surface, i.e. 

u ~ u thil e z . (9) 

In the surface integrals of eq. (|7) we start by considering 
the "front and back" surfaces perpendicular to the fins. Due 
to the translational symmetry along the fins and the vanish- 
ing of all gradients in this direction the integrals originating 
from the shear-stress contribution vanish and those orig- 
inating from the pressure contribution cancel when sum- 
ming over the front and back surfaces. The same applies to 
the surfaces parallel to the fins, since these faces constitute 
symmetry planes. Furthermore, since the flow is parallel to 
the top and bottom faces the integrals containing the pres- 
sure vanish, and we are left with the shear-integrals over 
these surfaces. From eqns. dSJ and © we find for the top 
surface 

r) I dA^UjdiUj - UjdiUj) = -rjUthrfA, (10) 

J top 
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where A is the area of the surface. The integral over the 
bottom surface only has contributions from the liquid-gas 
interface since the velocity vanishes on the fins 



T] I dA^UjdiUj—UjdiUj 

I bottom 



rr I dA y U z , 
liquid— gas 

(ID 

where we have used the fact that the reference flow only 
has components in z-direction and a vanishing shear rate at 
the liquid-gas interface. We thus arrive at 

Uth,i = —~a I dA y u z = p u (12) 

'// J liquid— gas V 

where the last equality comes from another application of 
the Lorentz reciprocal theorem taking a Couette flow with 
no-slip boundary condition everywhere on the bottom sur- 
face and shear rate 7 as "unhatted" flow, i.e. u = ■yye z , 
while retaining the Couette flow over longitudinal fins as 
reference flow. Also here only the integrals on the top and 
bottom surface containing the shear rate contribute, giving 
rj^PiA on the top surface and 777 j liquid _ gas dA y u z at the 
bottom. 

Equation ( fl~2l) is our result for the longitudinal thermo- 
capillary velocity in the Stokes limit. The corresponding 
longitudinal slip length is IToll 

2 ira 
Pi = In cos— , (13) 

7T Z 

where a = B/L is the liquid-gas fraction of the surface 
and Pi = p~i/L the dimensionless slip length in terms of 
half the fin spacing (c.f. figure [Tj. For later comparison, 
we invert relation ( PT2l) and define the dimensionless "lon- 
gitudinal thermocapillary slip coefficient" as 

T)/L 

PthA = u th,l Brr I n m\ - 

For realistic values of 10 wide fins with 40 fim spac- 
ing in between and an applied temperature gradient of 
(d z T) = — lOK/cm we obtain from equations ( [T3l and 
(fT2l . using the values of table J] 

P l = 0.7476, Pi = 18.7 fxm,u m = 2.9mm/s. (15) 

Note that this corresponds to a Reynolds number of approx- 
imately Re ~ 0.07, so the Stokes equation is still expected 
to be a good approximation. 

In the derivation above we assumed that the top surface 
of our control volume is located sufficiently high above the 
structured surface such that the thermocapillary flow field 
corresponds to a core plug flow. In fact, our numerical cal- 
culations confirm that this assumption is fulfilled reason- 
ably well already at a distance of 2L. 

Transverse fins - In the case of thermocapillary convec- 
tion perpendicular to the fins, the temperature gradient on 
the liquid-gas interface is not constant in the main flow di- 
rection and it is not possible to pull d z a out of the inte- 
gral, as done in equation ( TTTT ) . In fact, only in the case 



where the Marangoni stresses on the liquid-gas interface 
can be decomposed into a constant (leading to a main flow) 
and a part antisymmetric with respect to reflection at a fin 
center (contributing only to convective rolls but no main 
flow) can the Lorentz reciprocal theorem be applied as be- 
fore. Nevertheless, we expect the transverse slip coeffi- 
cient Pi = Pt/2 obtained for Couette flow over a finned 
geometry [10] to define an upper bound for the thermo- 
capillary slip velocity that can be obtained. In analogy to 
equation (fT4l we therefore define the "transverse thermo- 
capillary slip coefficient" 

•qa/L 

Pth,t = %,ts^ T y (16) 

In this definition we have incorporated the fact that 
the gradient driving the thermocapillary convection is 
(d z T) f ree = a -1 (d z T) since the temperature on the solid- 
liquid interfaces is approximately constant, owing to the 
assumed high thermal conductivity of the substrate. This 
seems to contradict the constant temperature gradient as- 
sumed in the bulk of the substrate. However, for fins much 
taller than wide, as usually applied for superhydrophobic 
surfaces, the diffusional character of the heat conduction 
equation assures that only the average temperature from the 
bottom prevails at the top of the fin. 

At this point a remark on the constant-temperature 
boundary condition on the solid-liquid interfaces is in or- 
der. This approximation will only be valid if the heat flux 
from the liquid is small enough. From dimensional con- 
siderations the wall heat flux in a fin scales as Jf in ~ 
pc p u t h,i{d z T)H / {\ — a). For an estimate we insert the 
values used to arrive at equation ( fT5l l together with a geom- 
etry height H = 200 /im and obtain as a measure for the 
temperature gradient in the fin Jfi n /ksi ~ 0.8K/cm <C 
{d z T), where the thermal conductivity of silicon ksi = 
148W/(mK) was used. Thus the temperature gradient in 
the fin is expected to be much smaller than the applied tem- 
perature gradient. The largest temperature gradient we will 
consider is (d z T) = — 60 K/cm, which from the above es- 
timate is seen to push our model to its limits. Nevertheless, 
this only gives a typical size of the temperature gradient in 
the fin, which will be mostly directed along its height and 
not along the solid-liquid interface. 

Numerical results - To test the analytical result (fT2l 
and to investigate the influence of inertia and non-uniform 
temperature gradients at the liquid-gas interface we have 
solved equations (O [6]) in the longitudinal case and equa- 
tions (O HI in the transverse case using a finite element 
discretization (Comsol mutiphysics). Due to symmetry 
a 2D calculation suffices in both cases and we restrict 
our attention to the liquid while assuming the tempera- 
ture of the substrate as given. The corresponding geome- 
tries are shown in figure |2 For definiteness we consider 
a fin spacing of 2L = 50 /xm and a geometry height of 
H = 200 fim while varying the liquid-gas interface area 
fraction a = B/L. The boundary conditions as well as 
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TABLE I: Boundary conditions and model parameters. The corresponding geometry is displayed in figure[2] 
boundary longitudinal transverse 

1, 1A, IB u x = u y =u z =0, T = T u y = u z = 0, T\ 1A = T , T\ 1B = T + 2L(d z T) 

2 T X y = -§frd X T, Ty Z = -§fr{d Z T),Uy=Q, dyT = Q Ty Z = - Sfr (d > Z T) , Uy = , dyT = 

3, 3A, W u x = 0, r xy = 0, d x u z = 0, d x T = u t \ 3A = Ui\ 3B , p\ 3A = p\sb, T yz = 0, 

r| 3 fl = T\ 3A + 2L{d z T), d z T\ 3B = d z T\ 3B 

4 Uy = 0, T X y = 0, dyU Z = 0, 9 y T = Uy = 0, Ty Z = 0, = 

77 = ImPas, p = 1000kg/m a , k = 0.6 W/(mK), c p = 4200 J/(kgK), J£ = -0.155 mN/(mK), T = 300 K 
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FIG. 2: (color online) (a) Simulation domains for calculating longitudinal (left) and transverse flow (right). The z-axis is always chosen 
along the main temperature gradient, the y-axis always normal to the structured surface. Boundary conditions and model parameters 
are specified in tableU (b) Longitudinal and (c) transverse thermocapillary slip coefficients as defined in eqs. d 1411 16b for different free 
surface fractions a — B/L and temperature gradients parametrized by (d z T) — — 10~-^ 60K/cm. The solid lines show the respective 
slip coefficients as calculated by Philip 1 10], i.e. eq. d!3b for the longitudinal case and half of this value for the transverse case. 



model parameters are given in table [IJ In particular, the 
wall temperature is prescribed, while the liquid-gas inter- 
faces are assumed to be adiabatic. 

Figure |2 shows longitudinal and transverse thermocapil- 
lary slip coefficients for different free surface fractions and 
temperature gradients. The solid lines show the respec- 
tive slip coefficients as calculated by Philip llfJl . In the 
longitudinal case slight deviations from the value obtained 
in the Stokes limit are observed at temperature gradients 
above (d z T) = -6 K/cm. This is expected, since for the 
largest temperature gradient of —60 K/cm the velocities are 
~ 25 mm/s for a = 0.9, corresponding to Re ~ 1. In the 
transverse case, some larger deviations appear. However, 
the data points stay below the limiting value obtained from 
the results of Philip even at low Re. Owing to non-constant 
temperature gradients along the liquid-gas interfaces, the 
flow velocity is reduced compared to the limiting value, 
an effect that shows up especially at larger Reynolds num- 
bers. Apparently, with the longitudinal arrangement signif- 
icantly larger flow velocities can be reached than with the 
transverse one. 

Conclusion - We have analytically derived a relation for 
the thermocapillary flow velocity along a finned superhy- 
drophobic surface. Even at moderate temperature gradients 
the resulting velocity values are large enough to suggest 
using this principle for microfluidic pumping. The pre- 
sented relationships between the slip length and the ther- 
mocapillary flow velocity give a very good approximation 



for the bulk fluid transport in the case of longitudinal fins 
and represent an upper limit to the flow velocity achiev- 
able with transverse fins. We expect that in the same spirit 
simple relationships for the thermocapillary flow velocity 
along alternative superhydrophobic surface structures can 
be found. 
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